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Abstract 



We establish universality of local eigenvalue correlations in unitary random matrix en- 
sembles -^\detM\ 2a e- ntlV ^dM near the origin of the spectrum. If V is even, and if the 

recurrence coefficients of the orthogonal polynomials associated with \x\ 2a e~ nVtyX ^ have a 
regular limiting behavior, then it is known from work of Akemann et a!., and Kanzieper and 
Freilikher that the local eigenvalue correlations have universal behavior described in terms 
of Bessel functions. We extend this to a much wider class of confining potentials V. Our 
(f) ' approach is based on the steepest descent method of Dcift and Zhou for the asymptotic 

analysis of Riemann-Hilbert problems. This method was used by Deift et al. to establish 
"p^ ' universality in the bulk of the spectrum. A main part of the present work is devoted to the 

|. analysis of a local Riemann-Hilbert problem near the origin. 

J: 

1 * 

1 Introduction 

In the present paper we consider the following unitary ensemble of random matrices, cf. [2] 

X i — \detM\ 2a e- ntTV{M UM, a > -1/2. (1.1) 

The matrices M are n x n Hermitian and dM is the associated flat Lebesgue measure on 
the space of n x n Hermitian matrices, and Z n is a normalizing constant (partition function). 
The confining potential V in (|1.1|) is a real valued function with enough increase at infinity, for 
example a polynomial of even degree with positive leading coefficient. Random matrix ensembles 
are important in many branches of mathematics and physics, see the recent survey paper 18 . 
The specific ensemble Ijl-ljl is relevant in three-dimensional quantum chromo dynamics 39 . 

The ensemble (jl.lj) induces a probability density function on the n eigenvalues x\, . . . , x n of 
M, given by 



1 n 
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Rn,m(yi,---,ym) = 7 — —^7 / •••/ P^ (Vl, ■ ■ ■ > 2/m, ^m+l , • • • , x n )dx m+ i . . . dx r . 



where Z n is the normalizing constant, and where w n is the following varying weight on the real 
line, 

Wn (x) = \x\ 2a e~ nV( - x \ for i G M. (1.2) 

Of particular interest are the local correlations between the eigenvalues of the ensemble (11.1(1 
of random matrices, when their size tends to infinity, at the origin of the spectrum, see |2"l l22llH2"] . 

The correlations between eigenvalues can be expressed in terms of the orthonormal polyno- 
mials Pk,n{x) = 'Jk,nX k + • • • with 7fc >ra > with respect to w n , that is 

Pk,n( x )PjA x )\ x \ 2ae ~ nV{x)dx = S jk- 

Namely, for 1 < m < n — 1, the m-point correlation function 

n\ 

(n — m)! 

V V ' 

n—m 

satisfies, by a well-known computation of Gaudin and Mehta |29| . 

l<i,j'<m) 

where 

n— 1 

K n (x,y) = \J w n (x)\/w n (y) ^Pj,n(x)p j>n (y) 

I T\ I 7-^7n-l,n Pn,n( x )Pn-l,n{y) ~ Pn-l,n( x )Pn,n(y) n Q , 

= y/ w n{x)Vw n (y) ! ; '■ '■ ; , (1.3) 

ln,n X-y 

which gives the connection with orthogonal polynomials. The second equality in (|1.3|) follows 
from the Christoffel-Darboux formula |35| . 

Akemann et al. j^l showed that the local eigenvalue correlations at the origin of the spectrum 
have a universal behavior, described in terms of the following Bessel kernel 

T a {u,v) = nV^V^^ ^ ^ , 1-4 

2{u — v) 

where J a ±i denotes the usual Bessel function of order a ± ^. In [2. it was assumed that the 
parameter a is a non-negative integer, that the potential V is even, and that the coefficients 
Cfc n in the recurrence relation 



x Pk,n{x) — Ck+l,nPk+l,n + c k,nPk- 



satisfied by the orthonormal polynomials have a limiting behavior in the sense that the limit 
Ck t n exists whenever k, n — > oo such that k/n — > t for some t > 0. The restriction that a is 
a non-negative integer was removed by Kanzieper and Freilikher [22, but they still required 
the assumption that V is even and that the recurrence coefficients have a limiting behavior. In 
fact, their method of proof (which they call Shohat's method) relies heavily on these recurrence 
coefficients. 
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It is the goal of this paper to establish the universality of the Bessel kernel (|1.4|> at the origin 
of the spectrum without any assumption on the recurrence coefficients. We can also allow V to 
be quite arbitrarily. We assume the following 

V : R -» R is real analytic, (1.5) 
V(x) 

lim ——\^— = +oo, (1.6) 

|x|-K» log(x 2 + 1) 

r/>(0) > 0, (1.7) 

where t/j is the density of the equilibrium measure in the presence of the external field V, |11[ I34]. 

Let us explain the condition (|1.7|) . Denote the space of all probability measures on R by 
Mi(R), and consider the following minimization problem 



inf / / log -dfi(s)d^t) + / V(t)dn(t) . (1.8) 

/Lt6Mi(R) \J J \S - t\ J ) 

Under the assumptions Q1.5|) and (|1.6f) it is known that the infimum is achieved EH uniquely 
at the equilibrium measure \iy € Mi(R) for V. The measure //y has compact support, and 
since V is real analytic, it is supported on a finite union of intervals. In addition it is absolutely 
continuous with respect to the Lebesgue measure, i.e. 

d(iy{x) = ip(x)dx, 

and tp is real analytic on the interior of the support of fj,y, see The importance of the 

equilibrium measure lies in the fact that ip is the limiting (as n — > oo) mean eigenvalue density 
of the matrix ensemble (|1.1[) , cf . 112j . The condition (|1.7f) then says that the mean eigenvalue 
density tp should be strictly positive there. If the origin belongs to the interior of the support 
of \xy but the mean eigenvalue density vanishes there, then the potential is called multicritical, 
see |31EJI^- This case will not be treated in this paper. 

The regular behavior of the recurrence coefficients assumed in [2~ll22| is probably satisfied if V 
is even and if the support of fj,y consists of one single interval. Note that we make no assumptions 
on the nature of the support of fiy. It can consist of any (finite) number of intervals. 

Our main result is the following. 

Theorem 1.1 Assume that the conditions ()1.5[) - (|1.7|) are satisfied. Letw n be the varying weight 
(jl.2|) . let K n be the kernel (|1.H[) associated with w n , and let i/j be the density of the equilibrium 
measure for V. Then, for u,v £ (0, oo), 

-K n ( =J° a (u,v)+0[ — V asn^oo, (1.9) 



m/>(0) n \ni/)(0y mp(0)J ay ' ' \ n 

where J° is the Bessel kernel given by (|1.4|) . The error term in ()1.9|) is uniform for u, v in 
bounded subsets of (0, oo). 

Other types of universal correlations have been established in the bulk [HI El I22I OH] , at the 
soft edge of the spectrum El E2I EUl EE] , and at the hard edge [HI EHl EH EH The universality 
at the hard edge is also described in terms of a Bessel kernel, which we have denoted in |28j by 
J Q , namely 

2{u — v) 

To distinguish with this Bessel kernel, we use J° to denote the Bessel kernel (jl.4|) relevant at 
the origin of the spectrum. 
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Remark 1.2 The universality is restricted to u,v > 0. It can be extended to arbitrary 
real u and v in the following way. For «,ti£l, we have that 

M-^- a ^ K "(^-^k)) =»-""-X(»,»)+0(i), asn^oo, (1.10) 

and the error term holds uniformly for u, v in compact subsets of R. We will restrict ourselves 
to proving but the same methods allow us to establish (jl.K)j) . 

Our proof of Theorem 11.11 is based on the characterization of the orthogonal polynomials 
via a Riemann-Hilbert problem (RH problem) for 2x2 matrix valued functions, due to Fokas, 
Its and Kitaev and on an application of the steepest descent method of Deift and Zhou 
|14j . See [HI HI] f° r an introduction. The Riemann-Hilbert approach gives asymptotics for the 
orthogonal polynomials in all regions of the complex plane, and it has been applied before on 
orthogonal polynomials by a number of authors, see for example [51 ^1 H31 123 EZ1 EE] • Bleher 
and Its |Hj and Deift et al ^2] were the first to apply Riemann-Hilbert problems to universality 
results in random matrix theory. Later developments include (SJ 1191 128j . 

In this paper we use many of the ideas of [12] • That paper deals with the varying weights 
e -nV(x) w [^Yl y satisfying Q1.5JI and Q1.6JI . The steepest descent method for Riemann-Hilbert 
problems is used to establish universality of the sine kernel in the bulk of the spectrum for 
the associated unitary matrix ensembles. In our case the general scheme of the analysis is the 
same, and we refer to [1 121 113j for some of the details and motivations. The extra factor \x\ 2a in 
our weights \x\ 2a e~ nV<yX ^ gives rise to two important technical differences. The first difference 
lies in the construction of the so-called parametrix for the outside region. To compensate for 
the factor |x| 2a we need to construct a Szego function on multiple intervals associated to \x\ 2a . 
The second and most important difference lies in the fact that we have to do a local analysis 
near the origin. This is where the Bessel functions J ,i come in. The construction of the 
local parametrix near the origin is analogous to the construction of the parametrix near the 
algebraic singularities of the generalized Jacobi weight, recently done by one of us in [SB]- The 
local parametrix determines the asymptotics of the orthonormal polynomials near the origin, 
and thus also governs the universality at the origin of the spectrum. 

The rest of the paper is organized as follows. In Section [2.11 we characterize the orthogonal 
polynomials via a RH problem, due to Fokas, Its and Kitaev [Ej. Via a series of transformations, 
we perform the asymptotic analysis of the RH problem as in |121 113j . The first transformation 
will be done in Section 12.21 the second transformation in Section [31 Next, we construct the 
parametrices for the outside region and near the origin in Section 0] and Section [SI respectively. 
The final transformation will be done in Section [13 Then we have all the ingredients to prove 
Theorem II .11 in Section Here we use some techniques from |28j . 

2 Associated RH problem and first transformation Y \— > T 

In this section we will characterize the orthonormal polynomials pk, n with respect to the weight 
(|1.2|) as a solution of a RH problem for a 2 x 2 matrix valued function Y(z) = Y(z;n,w), due 
to Fokas, Its and Kitaev [Ej) an d do the first transformation in the asymptotic analysis of this 
RH problem. 

2.1 Associated RH problem 

We seek a 2 x 2 matrix valued function Y that satisfies the following RH problem. 
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RH problem for Y: 

(a) Y : C \ R -» C 2x2 is analytic. 

(b) Y possesses continuous boundary values for x G R \ {0} denoted by Y + (x) and Y-(x), 
where Y + {x) and denote the limiting values of Y(z') as z' approaches x from above 
and below, respectively, and 

y + (x)=y_(x)^ |X| \ j, for*GR\{0}. (2.1) 

(c) Y(z) has the following asymptotic behavior at infinity: 
Y{z) = (l + o(-X\( Z l °V asar-oo. (2.2) 



*" 



(d) Y(z) has the following behavior near z = 0: 

r(*) = 



i |z' 2<: 
l u 



0\, !! 2 a)> if«<0, 



1 1 



O I J I ) , if a > 0, 



(2.3) 



as z -► 0, 2 G C \ R. 



Compared with the case of no singularity at the origin, see ^2]> we now have an extra 
condition l|2.3j) near the origin. This condition is used to control the behavior near the origin, 
see also |271I3B|. 



1 \z\ 2a 
\z\ 2a 



Remark 2.1 The O-terms in (j2.3|) are to be taken entrywise. So for example Y{z) = O ( ^ 
means that Yii(z) = O(l), r 12 («) = 0(M 2a ), etc. 

The unique solution of the RH problem is given by 

y(.) = 2m ^,„J-i (2.4) 

V-27r»7n-l,«Pn-l,n W "7n-l,n /i E^^l dx ) 

where pfc jn is the fcth degree orthonormal polynomial with respect to the varying weight w n , and 
where ^k,n is the leading coefficient of the orthonormal polynomial Pkn- The solution (|2.4|) is 
due to Fokas, Its and Kitaev pl|, see also [HI [T2l [TH1 - See [SUEZ] for the condition (TOl) . 

Note that (|2.4j) contains the orthonormal polynomials of degrees n—1 and n. By Q1.3JI it is 
then possible to write K n in terms of the first column of Y. So in order to prove Theorem ll.il 
an asymptotic analysis of the RH problem for Y is necessary. Via a series of transformations 
Y — > T ^ S ^ R we want to obtain a RH problem for R which is normalized at infinity (i.e., 
R(z) —* / as z —* oo), and with jumps uniformly close to the identity matrix, as n —* oo. Then i? 
is also uniformly close to the identity matrix, as n —* oo. Unfolding the series of transformations, 
we obtain the asymptotics of y. In particular, we need the asymptotic behavior of Y near the 
origin, which follows from the parametrix near the origin. 
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2.2 First transformation Y — ► T 

We first need some properties of the equilibrium measure \iy for V. Its support is a finite union 
of disjoint intervals, say [Jj^i[bj-i, a j\- S° the support consists of N + 1 intervals and we refer 
to these as the bands. The complementary N intervals (a,j,bj) are the gaps. Following ^2j, we 
define 

N+l 

J= U ( 6 i-i' a i) 

so that J is the interior of the support. The density ip of fly has the form |12| 

^(x) = — R^fxWx), for x G J, (2.5) 

2lTl 



where 



N+l 

R(z)=H(z-b J - 1 )(z-a j ), (2.6) 

3=1 

and where ft, is real analytic on R. In this paper we use 72 1//2 to denote the branch of \fH which 
behaves like z N+1 as z —> oo and which is defined and analytic on C \ J. In (|2.5|) we have that 
R 1 ^ 2 denotes the boundary value of R l l 2 on J from above. There exists an explicit expression 
for h in terms of V, see jllj . but we will not need that here. 

The equilibrium measure minimizes the weighted energy (|1.8j) . The associated Euler-Lagrange 
variational conditions state that there exists a constant I G R such that 

2 / log Is - s\i(j(s)ds - V{x) = £, for x G J, (2.7) 



2 / log |x - s\ip(s)ds - V(x) < £, for x G R \ J. (2.8) 

The external field is called regular if the inequality in ()2.8|) is strict for every x G R \ J, and 
if /i(x) 7^ for every x G J. Otherwise, V is called singular. The regular case holds generically 
|25| . In the singular case there are a finite number of singular points. Singular points in J are 
such that h vanishes there. Singular points in R \ J are such that equality holds in (|2.8|) , 

In order to do the first transformation, we introduce the so-called (7-function ^21 Section 3.2] 



g(z) = J \og{z - s)tp(s)ds, for z G C \ (-00, a N+1 ], (2.9) 

where ip(s)ds is the equilibrium measure for V . In p. 9)1 we take the principal branch of the 
logarithm, so that g is analytic on C \ (—00, azv+i]- 

We now give properties of g which are crucial in the following, |12| Section 3.2]. From the 
Euler-Lagrange conditions ()2.7() and Q2.8|) it follows that 

g+(x) +g-(x) -V(x) -£ = 0, for x G J, (2.10) 

g + (x)+g-(x)-V(x)-£<0, for x G R \ J. (2.11) 

A second crucial property is that 

ra-N+l 

g + (x) - g-(x) = 2m I d^ v (s), for x G (-00, a N+1 ), (2.12) 
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so that g+(x) — g~{x) is purely imaginary for all x £ R and constant in each of the gaps, namely 



5+0*0 -5-0*0 



' 2vri, 

2irif^ N+1 dfj, v (s) =: 2vri%, 
> 0, 



for x < bo, 

for x £ {cLj , bj), j = l 
for x > 0/v+i- 



(2.13) 



From (|2.13j) we see that £lj is the total ^y-mass of the N + l — j largest bands. These constants 
all belong to (0, 1). Note that Oj was defined with an extra factor 2tt in |12j . 

As in |12l Section 3.3], we define the matrix valued function T as 



T(z) = e-- (73 Y(z)e-' T3 e 



for z e c \ : 



(2.14) 



where 03 = ( o -i ) i s the Pauli matrix. Then T is the unique solution of the following equivalent 
RH problem. 

RH problem for T: 

(a) T : C \ R -» C 2x2 is analytic. 

(b) T satisfies the following jump relations on R: 



T+(x) = T_(x) 



-n(fl+(as)-p_(a;)) 




X 



2a 



n(g + (x)-g-{x)) ) ' 



for x G J\ {0}, 



(2.15) 



T+(x)=T_(x) 



g—2mnQj 



\ x \2a e n(g + (x)+g_(x)-V (x)- 
^2ninflj 



(2.16) 



for x e (aj,bj),j = 1...,N, 



(\ \ x \^ e n(g + (x)+g-{x)-V(x)-i)\ 

T + (x) = T-(x) I 1 1 , for x < b or x > o/v+l- ( 2 -17) 

(c) T(z) = J + 0(l/z), as z -» oo. 

(d) T(z) has the same behavior as Y(z) as z — > 0, given by (|2.3I) . 

3 Second transformation T — > 5 

In this section we transform the oscillatory diagonal entries of the jump matrix in (|2.15[) into 
exponentially decaying off-diagonal entries. This lies at the heart of the steepest descent method 
for RH problems of Deift and Zhou [Hj, and this step is often referred to as the opening of the 
lens. 

For every z£C\M lying in the region of analyticity of h, we define 

<j>{z) = - / RV 2 {8)h(8)d8 
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where the path of integration does not cross the real axis. Since R 1 ^ 2 (s)h(s)ds = for 
every k = 1,...,N, (this follows easily from the formulas in |12l Sections 3.1 and 3.2]), and 
j-ajv+i R 1 ^ 2 (s)h(s)ds = 2iriQ,j, we find that for every j, 

/a,j 
R 1/2 {s)h(s)ds + 27rin jl iflmz>0, (3.1) 

R 1/2 {s)h(s)ds-27rin jl iflmz<0. (3.2) 



Note that in |12j a function G is defined which is analytic through the bands. We found it more 
convenient to have a function with branch cuts along the bands, see also [H]. The functions G 
and <p also differ by a factor ±2. 

The point of the function (j) is that 4> + and 4>- are purely imaginary on the bands, and that 

20 + = -2</>_= 5+ - 5 _. (3.3) 

This means that 2cp and — 2<p provide analytic extensions of g+ — g_ into the upper half-plane 
and lower half-plane, respectively. We also have that for z in a neighborhood of a regular point 
x £ J, (see ^1 Section 3.3] for details) that 

Re<j)(z)>0, iflmz/O, (3.4) 

We will now discuss the opening of the lens in the regular case. In the singular case we need 
to modify the opening of the lens somewhat, since we have to take into account the singular 
points. We do not open the lens around singular points that belong to J, see |12l Section 4] for 
details. 

For V regular, there is a suitable neighborhood U of J such that the inequality in Q3.4JI holds 
for every z S U. The opening of the lens is based on the factorization of the jump matrix ()2.15|) 
into the following product of three matrices, see also (|3.3j) . 

' e - n (g+( x )-g-( x )) \x\ 2a \ f e ~ 2n( t>+{ x ) \x\ 2a 

e n(g+(x)~g-(x)) J ~ I q e -2n0_(z) 

1 0\ ( \x\ 2a \ ( I U\ 

(3.5) 



^ x ^-2a e -2ncf)-(x) -^J ^ j^-j — 2ct q J \^ x \~2a e -2ncj>+{x) ^ 

As in |38) we take an analytic contination of the factor \x\ 2a by defining 

f (-z) 2a , ifRez<0, 
io(z) = \ \ (3.6) 
\ z 2a , if Rez > 0, 

with principal branches of powers. In contrast to the situation in [2], here we have to open the 
lens also going through the origin, cf. j^Hj- This follows from the fact that \x\ 2a does not have 
an analytic continuation to a full neighborhood of the origin. 

We thus transform the RH problem for T into a RH problem for S with jumps on the oriented 
contour S, shown in Figure ^ The precise form of the lens is not yet defined, but it will be 
contained in U. 

Define the piecewise analytic matrix valued function S as 

T(z), for z outside the lens, 

.S'(-) { T ( z ) (_ w ( z )-i e - 2n< «^ i) ' for z in the upper parts of the lens ' (3.7) 
T{z) ( ^ _2n<t)(z) > ^ or z m * ne l° wer parts of the lens. 

Then, S is the unique solution of the following equivalent RH problem. In (|3.8|) . C+ and C_ are 
used to denote the upper half-plane {Imz > 0} and the lower half-plane {Im z < 0}, respectively. 
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Figure 1: Part of the contour S. 



RH problem for S: 

(a) S : C \ S -» C 2x2 is analytic. 

(b) S 1 satisfies the following jump relations on £: 

S + (^)=fir_( z )^ ) _ 1 ^ w JV forzGSnC ± , (3.8) 



,S'_(.r) ,S'_(.r)( ° ^l 2 ^) , (or.rG.A {()}. (3.9) 



■ „—2irinQj \ q ,\2a p n(g+(x)+g-(x)—V(x)—£)' s 

S + (x) = S^(x)[ e Q |X| g2mnCj j. (3.10) 



for a; £ (aj,bj),j = 1. . .N, 



S+(x) = S-(x) ( j J, for x < 6o or a; > ajv+i- (3.11) 

(c) S(z) = / + 0(1/2), as z -» 00. 

(d) For a < 0, the matrix function S(z) has the following behavior as z — > 0: 

= O ' asz^0,zGC\S. (3.12) 



For a > 0, the matrix function 5"(z) has the following behavior as z — > 0: 

5(z) = < 



O ( j jj, as z — > from outside the lens, 



(3.13) 



() [ I ^ -2a ]_ ) > as z ^ from inside the lens. 

(e) S remains bounded near each of the endpoints Oj, bj. 

By (|3.4|) the factor e _2n9i ( 2 ) in (|3.8() is exponentially decaying for z 6 S n C± as n — ► 00. 
This implies that the jump matrix for S converges exponentially fast to the identity matrix as 
n — ► 00, on the lips of the lens. Since V is regular, we have the strict inequality 

g+(x) +g-(x) -V(x) -£ < 0, forxGM\J, (3.14) 

so that the factor e n ^ 9+ ^ +9 ^^~ v ^ x ^^ in ()3.1(Jj) and ()3.11j) is also exponentially decaying as 
n — > 00. 
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4 Parametrix for the outside region 

From the discussion at the end of the previous section we expect that the leading order asymp- 
totics are determined by the solution of the following RH problem. 

RH problem for p(°°) : 

(a) P(°°) : C\ [b ,a N+1 ] -► C 2x2 is analytic. 

(b) p(°°) satisfies the following jump relations: 

Pt\x) = Pi°°\x) (_ |a J_ te |X f), forxGJ\{0}, (4.1) 

^^"'W 6 e27r L,J, forxG(a„6,),i = l...iV, (4.2) 

(c) P(°°)(z) =7 + (1/z), as z ->■ oo. 

The solution of this RH problem is referred to as the parametrix for the outside region, and will 
be constructed using the so-called Szego function on the union of disjoint intervals J, associated 
to \x\ 2a . The importance of the Szego function is that it transforms this RH problem into a RH 
problem with jump matrix ( _?i o) on J. 

4.1 The Szego function 

We seek a scalar function D : C \ [bo, o-n+i] — ► C that solves the following RH problem. 

RH problem for D 

(a) D is non-zero and analytic on C \ [bo, ajv+i] 

(b) D satisfies the following jump relations: 

D + {x)D_(x) = \x\ 2a , forx£j\{0}, (4.3) 
D+(x) = e 27r ^D_(x), for x G ( aj , bj), j = 1, . . . N, (4.4) 

for certain unknown constants £i, . . . , £jv e ^- The selection of £i, . . . , £jv is part of the 
problem. We should choose them such that it is possible to construct D. 

(c) D and D^ 1 remain bounded near the endpoints ai,bj of J, and 

Doo := lim D(z) (4.5) 

z^oo 

exists and is non-zero. 

We seek D in the form D(z) = exp$(z). Then the problem is reduced to constructing a 
scalar function <I>, analytic on C \ [feojO/v+i], remaining bounded near the endpoints a^, bj of J 
and at infinity, and having the following jumps 

+ 3>_(x) = 2a log |x|, forxGJ\{0}, (4.6) 

$ + (a;) = $_(x) + 27ri&, forx€(oj,6j), j = l,...,JV. (4.7) 
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We can easily check, using Cauchy's formula, the Sokhotskii-Plemelj formula [201, an d the 
fact that R l _l 2 {x) = -R l / 2 (x) for x G J, see (ESJ), that $ defined by 



K) 1 ^ fl i/2 (:B) x-; 



TV 



, V •/„, flV*(s)s- 



(4.8) 



satisfies the jump conditions Ij4.fi|) and (|4.7|) . We note that $ is analytic on C \ [&o,o/v+i] and 
remains bounded near the endpoints cn,bj of J. We use the freedom we have in choosing the 
constants £i, . . . ,£jv to ensure that 3> remains bounded at infinity. Since R 1 ^ 2 {z) behaves like 
oo, and since 



N+l 

z ~ as z 



N-l u 
x K 



x — z 



E 

fc=0 



fc+1 



+ o 



1 



JV+l / ' 



as z —>■ oo, 



we have to choose £i, . . . , £jv such that the ./V conditions 
1 f 2a log I x I k , ^ x fc dx 



0, 



,iV-L 



(4.9) 



are satisfied. Note that (|4.9j) represents a system of N linear equations with coefficient matrix 

bi dx r&2 dx rb N dx 

Jao Ri72iVl ' ' ' Ja N W^Jxj 

fbN xdx 
Jc 



A 



(r 

Jai 

fbi xdx rOj 

Jai Rl/2(x) Ja 2 R}/l{x) 



R 1 / 2 (x) Ja 2 WT^x) 
l xdx f&2 xdx 



>a N 



(4.10) 



rbi x N ~ 1 dx r&2 x^" 1 ^ rbjv x^" 1 ^ . 

\Jai R 1 / 2 (x) J a,2 R 1 / 2 (x) " ' Ja N flV2(x)/ 

By the multilinearity of the determinant, we have 

/I 1 ••• 



det A 



Ul 



det 



ttjv 



1 \ 



N-l _JV-1 



dXN 



RV 2 ( X1 ) R l l 2 (x N ) 



d,XN 



x x ) R}/\ XN y 



(4.11) 



Since the gaps (aj, bj) are disjoint, we have Xj < x^ for j < k, so that Ylj^fak ~ x j) > 0- Using 
the fact that R 1 ! 2 does not change sign on each of the gaps (aj,bj), it then follows that the 
integrand in (|4.11j) has a constant sign in the region of integration, so that det A ^ and A is 
invertible. 

We then define £i, . . . , £jv as follows 



/6\ 
VW 



/ J_ r 2a log |x| , \ 
/ 2^ J J ^1/2^) aX \ 



1 r 2a 

2tt?. J J R i 



2a log \ x\ 



xdx 



1 f 2a log 1 1 1 JV-l , 
\2™Jj ^/2 (x) a 



(4.12) 
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Figure 2: The oriented contour 75. 



Note that R 1 / 2 is real on each of the gaps, so that by 1)4.10(1 all entries of A are real. This implies, 

1 /2 

from (|4.12f) and the fact that R / is purely imaginary on J, that the constants £1 , ■ ■ ■ , £jv are 
real. 

We proved the following 

Theorem 4.1 The scalar function D{z) = exp$(z), where <1> is given by (|4.8[) . and tae con- 
stants £1, • • • , £jv % <|4.12|) . solves the RH problem for D. 

For later use we state the following lemma. 

Lemma 4.2 We have that z~ a D(z) and z a D(z)~ 1 remain bounded near the origin. 

Proof. For definiteness, suppose that the origin lies on the band (bj, a^+i) with j E {0, . . . , N}. 
Since D(z) = exp &(z), with $ given by ()4.8() . it is sufficient to prove that 

1 P+ 1 2alog|x| dx a log z 

+ F(z), (4.13) 



2™ 4 R^ix)*-* R 1/2 {z) 

with .F analytic near the origin. Fix z near the origin with Im z 7^ 0, and let 75 be the oriented 
contour shown in Figure O with 5 > small. Cauchy's formula implies 

1 /" logC d( logz 
2^ 4 i2V2(C)C-^ = W\z)~ 

1/2 1/2 

Letting 5 — ► 0, we then have, since i2 + (x) = — i?_ (x) for x E (6j,a J+ i), 
~, . 1 f° log Ixl + in dx 1 loglxl dx 

F i z ) + ^~ —jtt—z — : + 



27ri -'bj R]/ 2 (x) x ~ z 2 ™Jo R 1 l 2 {x) x ~ z 



2 ™ Ja f+1 R\ 



g\x\ dx 1 f 3 log|x|— iir dx logz 



Miiiffi) 1 " 2 2 ™^° i?+ /2 (x) z-z fi 1 / 2 ^)' 

with F analytic near the origin. Hence 

1 /' a J+ 1 log|x| dx 1 logz 1 ~ 
2^4 i?f(x)^-^ = 2i^)"2 W ' 

so that 1)4.13(1 holds with F(z) = —aF(z), which proves the lemma. □ 
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4.2 Construction of P(°°) 

We now use the Szego function D from the previous subsection, to transform the RH problem 
for p(°°) into a RH problem with jump matrix ( _?i J) on J- We seek p(°°) in the form, cf. 

PW{z) = D2P (co) (z)D(z)-^, for z e C\[b ,a N+1 }. (4.14) 

Then, by (|4.1|) - (|4.4jl the problem is reduced to constructing a solution of the following RH 
problem. 

RH problem for p(°°) 

(a) P(°°) : C \ [b , a N+1 ] -► C 2x2 is analytic. 

(b) p(°°) satisfies the following jump relations: 

P^°\x) = P^\x) f ^ J J , for x G J \ {0}, (4.15) 



p} oo ) (x ) = J p(°°)( 2; )( e ; c ^ e2 ™4_ 2 ^), (4.i6) 



:inQj 

for x G (dj,bj), j = 1, . . . ,N. 



(c) p(°°)(z) = / + 0(l/z), as z -»■ oo. 

This corresponds to the RH problem |121 (4.24)-(4.26)], which has been solved there using 
Riemann theta functions. Note that, in contrast to the RH problem j!21 (4.24)— (4.26)], the 
jump matrix in Q4.16JI contains extra factors exp(±27ri£j) in the diagonal entries, which come 
from the Szego function D. However, this does not create any problems. 

In order to formulate the solution of the RH problem for p(°°) we need to introduce some 
additional notations. Here we closely follow ^2], see also |1U| . 

Let J = R \ J be the complement of J, and ao = ojv+l- Letting the point oo lie on the 
interval (ao,6o)> J can be displayed as a union of intervals on the Riemann sphere. Let X be 
the two-sheeted Riemann surface of genus N associated to \J R(z), obtained by gluing together 
two copies of the slit plane C \ J along J. We draw cycles Aj winding once, in the negative 
direction, around the slit (a,-, bj) in the first sheet, and cycles Bj starting from a point on the slit 
(dj,bj) going on the first sheet trough a point on the slit (ao,bo), and returning on the second 
sheet to the original point, as indicated in FigureEl The cycles {Ai, Bj}i<i,j<N form a canonical 
homology basis for X, see |15| . 

Let uj = (u>i, . . . ,ujf) be the basis of holomorphic one-forms on X dual to the canonical 
homology basis, that is 



Ia 



wi = <5y, l<i,j<N. (4.17) 
The associated Riemann matrix of B periods, denoted by r and with entries 

Tij= [ Ui, l<i,j<N, (4.18) 
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Figure 3: The canonical homology basis {A{, -Bj}x<ij<Ar for X. The full lines denote paths on 
the first sheet, while the dotted lines denote paths on the second sheet. 

is symmetric with positive definite imaginary part, see The associated Riemann theta 

function is defined by 



9{z) = exp 27ri ^(m, z) + — (m, rm)j , 



zeC ! \ (4.19) 

where (•, •) is the real scalar product, which defines an analytic function on C . The Riemann 
theta function has the periodicity properties ^3] with respect to the lattice 7L N + tZ n 

8(z + ej ) = 6(z), 9(z ± Tj ) = e T27Ti ^- niT ^8(z), (4.20) 

where z = (z%, . . . , z^) and e,- is the jth unit vector in with 1 on the jth entry and zeros 
elsewhere, and where Tj is the jth column vector of r. 
Define the scalar function 



TT ( z ~ b A ( z ~ b o \ 
iJl \z-a>i) \z- a N +i) 



1/4 



(4.21) 



which is analytic on C \ J, with 7(2) ~ 1 as z —* 00, z G C+. It is known |12l Lemma 4.1] that 
7 has the following properties: 

7 + 7 _1 possesses N roots {zj ^}jL\ with zj ^ on the — side of (aj, bj), 

7 — 7 _1 possesses ./V roots {z^}^ =l with on the + side of (aj, bj). 

Fix the base point for the Riemann surface X to be a^+i = ao ; let K be the associated 
vector of Riemann constants ^5] , and define the multivalued function 

u(z) = I uj. (4.22) 

J ajv+i 



Here, we take the integral along any path from a^+i to z on the first sheet. Since the integral is 
taken on the first sheet, u{z) is uniquely defined in C N /7* N because of (|4.17|) . Let d be defined 

as 



N „ z {~) 



d = -K-^2 (4.23) 

j = l J a N+l 

where again the integrals are taken on the first sheet. 
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We now have introduced the necessary ingredients to formulate the solution of the RH 
problem for p(°°). Together with Q4.14JI this gives the parametrix p(°°) for the outside region. 
The solution of the RH problem for p(°°) is given by, see |12| Lemma 4.3], 

P(°°)r^ Hintr ( d{u+(co)+d) e(u+(oo)+d) 

r \t) — U(u + (oo)-nSi+C+(i)' ^-u+tooJ-nSJ+C-! 



for z G C+, and 

p(°°)(z) = diag 



Py_-_1 0(u(z)~nn+£+d) 7-7- 1 fl(-n(z)-nfi+g+rf) ' 
2 0(u(.z)+<f) -2i 0(-u(z)+d) 

7-7- 1 e(w(.z)-nn+g-d) 7+7- 1 e(-u(z)-n^+g-rf) 
2i 0(u(z)-d) 2 6»(«(z)+d) 



6(u + (oo)+d) 6(u + (oo)+d) 
e(u + (oo)~nQ+(+d) ' 0(-u + (oo)-nf2+§-(2) 



-2i 9(-u(z)+d) 2 6>(u(z)+d) 

7+7- 1 6(-u(z)-nQ+£-d) 7-7- 1 6>(u(z)-nn+g-d) 
2 0(u(z)+d) 2i 0(u(z)-(i) 



(4.24) 



(4.25) 



for z G C_. Here, = (f2i, . . . ,J1jv) and £ = (£1, . . . , £jv) 



Remark 4.3 In contrast to |12j we have an extra term £ in the Riemann theta functions. This 
comes from the slightly different jump matrix in (|4.2jl due to the Szego function, as noted before. 
If £ G the factors e ±27r4 ^ in (|4.1fiJI disappear and the RH problem for p(°°) is exactly the 
same as the RH problem |12l (4.24)-(4.26)]. Since the Riemann theta functions possess the 
periodicity properties (|4.2()jl . the term £ in (|4.24j) and (|4.25[) disappears in this case. This is in 
agreement with |12l Lemma 4.3]. 

For later use, we need p(°°) to be invertible. In |12| Section 4.2] it has been shown that 
det P(°°) = 1, so that by gUP 

detP (oo) = l. (4.26) 

5 Parametrix near the origin 

In this section we construct the parametrix near the origin. As noted in the introduction, it is 
similar to the construction of the parametrix near the algebraic singularities of the generalized 
Jacobi weight [HE]) and we skip some details and motivations. 

We surround the origin by a disk Us with radius 5 > 0. We assume that 5 is small, so that 
in any case, we have that [—5, 5] C J. We seek a matrix valued function P that satisfies the 
following RH problem. 

RH problem for P: 

(a) P(z) is defined and analytic for z G U§ \ £ for some 5o > S. 

(b) On £ n Us, P satisfies the same jump relations as S, that is, 

p + (z) = P-(z)( u{z) J e - 2n4){z) °), forzGSn([/ 5 nC ± ), (5.1) 
P + (x) = P_(x)(_, ° for*G(-M)\{0}. (5.2) 
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(c) On dUs we have, as n — > oo 

P(s) (V (oo) ) _1 (z) = / + O (^) , uniformly for zedU s \ S. 

(d) For a < 0, the matrix function P(z) has the following behavior as z — > 0: 



(5.3) 



P(Z) = °( V 1 |z| 2 -J' as ^°- ( 5 - 4 ) 
For a > 0, the matrix function P(z) has the following behavior as z — > 0: 



P(z) = 



Of | | ) . as ; — () from outside the lens. 

M-2« 1^ 

^ I -2a i ) ' as z ^ from inside the lens. 



(5.5) 



We construct P as follows. First, we focus on conditions (a), (b) and (d). We transform 
the RH problem for P into a RH problem for P^ with constant jump matrices, and solve the 
latter RH problem explicitly. Afterwards, we also consider the matching condition (c) of the 
RH problem. 

We start with the following map / defined on a neighborhood of the origin 

f i0(z)-»0 + (O), iflmz>0, 
J\ z ) = \ (5.6) 
[ -i<M0), iflmz<0. 

Since = —</>_, we have that / is analytic for z in a neighborhood of the origin. An easy 
calculation, based on the fact that 2cj) + {x) = g+(x) — g~(x) and on ()2.9j) and (|5.6jl . shows that 



/(x) = vr / ^(s)ds, for x G (-5,5), (5.7) 
Jo 

which implies that /'(0) = 7r^(0) > 0. So, the behavior of / near the origin is given by 

f(z) = irip(0)z + O (z 2 ) , as z 0. (5.8) 

So if we choose 5 > sufficiently small, C = f( z ) 1S a conformal mapping on onto a convex 
neighborhood of in the complex £-plane. We also note that f(x) is real and positive (negative) 
for x £ Us positive (negative), which follows from (|5.7|) . 
Let Tj, j = 1, . . . , 8 be the infinite ray 

={(SC |argC ={j-i)\}- 

These rays divide the (-plane into eight sectors I- VIII as shown in Figure H We define the 
contours Ej, j = 1, 2, . . . , 8 as the preimages under the mapping ( = f(z) of the part of the 
corresponding rays Fj in f(Us), see FigureEl 

We have some freedom in the selection of the contour S. We now specify that we open the 
lens in such a way that 

znu s = J Zj. 

7=1,2,4,5,6,8 

As a consequence we have that / maps X to part of the union of rays (J ■ Tj . 
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In order to transform to constant jumps we use a piecewise analytic function W corresponding 
to the analytic continuation of |x| 2a . For z £ Us, we define 

if./2<|a rg/W |<., 
\ (-*)«, ifO<|«rg/W|<7r/2, 

with principal branches of powers. Then W is defined and analytic in Us \ (£1 U S3 U £5 U £7). 
We seek P in the form 

P{z) = E n {z)P^{z)W{zy a,i e- n ^ z)cTi . (5.10) 

Here the matrix valued function E n is analytic in a neighborhood of Us, and E n will be deter- 
mined below so that the matching condition (c) of the RH problem for P is satisfied. Similar 
considerations as in |38| show that P^ should satisfy the following RH problem, with jumps on 
the system of contours Ui=i oriented as in the left part of Figure El In (j5,ll|) - (|5.14j) . £° is 
used to denote £j without the origin. 

RH problem for P^: 

(a) (z) is defined and analytic for z G Us \ (£ U T) for some 6q > 5. 

(b) P^ satisfies the following jump relations on Us D (£ U T): 

forxe£?U£!?, (5.11) 

Jj, forzG^USg, (5.12) 
forzG£§U£?, (5.13) 
?), forze£°U£g. (5.14) 



P} 1} (x) 


= p(l) 




pf(.) 






P«(.) 


= 


[z)e maa3 , 


p«(.) 







(c) For a < 0, P^^(z) has the following behavior as z — > 0: 



p(1)(z) = °l|J uP' asz ^ a (5 - 15) 



For a > 0, PW(.z) has the following behavior as z — > 0: 
P«(z) = 



O ( ] L -a ) ' as z ^ from outside the lens, 



\ Z \— OL \%\~ a 

O ( _ a -a ] ) as z ^ from inside the lens. 



(5.16) 



Next we construct an explicit solution of the RH problem for P^\ This is based on a model 
RH problem for "if a in the C-plane, see [SHI - We denote by Tq, the contour U? = i oriented as 
shown in Figure |1] 
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RH problem for *& a : 

(a) ^ a : C \ -> C 2x2 is analytic. 

(b) ^ a satisfies the following jump relations on Tq,: 

*« 1+ (0 = ^.-(Of^ J), forCerxUTs, (5.17) 

* a>+ (0 = *«,-(C)( e -L, J). forCer 2 ur 6 , (5.18) 
* a ,+(C) = * a ,-(C)e™, forcer 3 ur 7 , (5.19) 
* a ,+(C) = vPa.-CofgaL ?)> forCGr 4 ur 8 . (5.20) 

(c) For a < the matrix function ^q(C) has the following behavior as £ — > 0: 

*«(C) = o(|^p asC-0. (5.21) 

For a > the matrix function ^f a {C) nas the following behavior as £ —* 0: 



*a(C) = < 



^(j^jl j^j-J- hsc -0 with , e II. III. VI. VII. 



n flCI~ Q Kr\ 

I vicr icrv 



(5.22) 



as C -> with C € I, IV, V, VIII. 



This RH problem was solved in |381 formulas (4.26)-(4.33)]. It is built out of the modified Bessel 
functions I„.i, K,i and out of the Hankel functions H , , , H^ 2 } , . For our purpose here, it 
suffices to know the explicit formula for ty a in sector I. There we have 

1 /i^Mo -^ (1) !(cA 

*a(C) = ^V^C 1/2 Z 2 m 2 e-( Q+ i)™3 ) forO<argC<f. (5.23) 

2 K4 (c) -<V c v 

Starting from (|5,23|) we can find the solution in the other sectors by following the jumps (|5.17|) - 
l|5.2U|) . See [HE! f° r explicit expressions. 
Now we define 

pW(z) = * a (nf(z)), (5.24) 

and P^ will solve the RH problem for P^\ This ends the construction of P^\ 

So far, we have proven that for every matrix valued function E n analytic in a neighborhood 
of Us, the matrix valued function P given by 

P{z) = E n (z)^ a (nf(z))W(z)-^e- n ^ z ^, (5.25) 

satisfies conditions (a), (b) and (d) of the RH problem for P. We now use the freedom we have 
in choosing E n to ensure that P, given by (|5.25|) . also satisfies the matching condition (c) of 
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the RH problem for P. To this end, we use the asymptotic behavior of VPq, at infinity, see |38| 
(4.43)— (4.46)]. Similar calculations as in [3S] show that we have to define E n as 

E n {z) = B( z )e^e>i Q j) , (5.26) 

where the matrix valued function I? is given by 

£(2) = P^(z)W(z) a3 e^ ani<73 , for z G / _1 (/U /I), (5.27) 

£(2) = p(°°)(z)^(z) a3 e-i a7ri,J3 , for z G f~ l (III U JV), (5.28) 





Ot7vi(T 3 

? 




i a7ri(T3 


CT3 ^ 


1 




-1 


0-3 f 


1 




-1 



E(z) = ^'(^(z) 13 ^ Ag-iTrio^ for z 6 /-i(VUVi), (5.29) 

= p(°°)(z)W(z) ff3 ^je^ ma(7i , for z G f-^VII U VIII). (5.30) 

Following the proof of [38( Proposition 4.5], we obtain that E is analytic in a full neighbor- 
hood of Ug. Here we need the fact that D{z)/W{z) and W(z)/D{z) remain bounded as z — > 0, 
which follows from (|5.9j) and Lemma 14.21 Then we see from (|5.26j) that E n is also analytic in a 
neighborhood of Ug. This completes the construction of the parametrix near the origin. 

Remark 5.1 Note that, in contrast to the case of the generalized Jacobi weight [SHI) here E 
depends on n. This follows from the fact that the parametrix for the outside region in our 
case depends on n. 



For later use we state, since E is analytic in Ug and from the explicit form of p(°°\ cf. |12j . 

A 

dz ' 



that E(z) and 4^E(z) are uniformly bounded for z G Ug, as n — > oo 



6 Third transformation S — > 

At each of the endpoints Oj, 6,- of J, we have to do a local analysis as well as at each of the 
singular points (if any). The endpoints and singular points are surrounded by small disks, say 
of radius 5, that do not overlap and that also do not overlap with the disk Ug around the region. 
Within each disk we construct a parametrix P which satisfies a local RH problem: 



RH problem for P near Xq where xq is an endpoint or a singular point: 

(a) P(z) is defined and analytic for z £ {\z — xo\ < 5} \ Y> for some 5q > 5. 

(b) P satisfies the same jump relations as S does on £ n {\z — xq\ < S}. 

(c) There is k > such that we have as n — > oo: 

P(z) 1 (z) = I + O (^\ , uniformly for \z - x \ = 5. (6.1) 

The local RH problem near the regular endpoints a,, bj of J is similar to the situation in |T2j. 
Here however, we have extra factors |x| ±2a and cj(z) _1 in the jump matrices. These factors can 
easily be removed via an appropriate transformation, and the local RH problem is then solved 
as m [T2l Section 4.3-Section 4.5] with the use of Airy functions. For our purpose, we do not 
need the explicit formulas for the parametrix near the endpoints. It suffices to know that P 
exists. For regular endpoints we can take k = 1 in l|6.1|) . 
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b H a h o a i+ i b 



Figure 6: Part of the contour T,r. The singular point Z\ corresponds to a point where h vanishes 
at the interior of J, the singular point z<i corresponds to a point where we obtain equality in 

m 



Near the singular points we can follow the analysis of |12) Section 5]. Here we do not 
construct an explicit parametrix out of special functions, only existence of the local RH problem 
is obtained. For singular points we have k < 1 in (|dlj) . see ^2]- So the singular points lead to 
error terms with decay slower than for the regular points (k = 1). However, this has no influence 
on the universality result at the origin (not even in the error term), since that only depends on 
the leading order asymptotics. 

We are now ready to do the final transformation. As noted before, we surround the endpoints 
aj, bj of J, the origin, and the singular points of the potential V by nonoverlapping small disks. 
Using the parametrix p(°°) for the outside region and the parametrix P defined inside each of 
the disks, we define the matrix valued function R as 

f S(z) (Pi 00 ))' 1 (z), for z outside the disks, , . 

R(z) = I V n _ ' V h ' (6.2) 

[ S(z)P (z), for z inside the disks. 

Remark 6.1 It is known that the inverses of the parametrices p(°°) and P exist, since all 
matrices have determinant one. For P(°°), see 1)4.26(1 . For P within the disks around the 
endpoints en, bj of J, as well as within the disks around the singular points of V we refer to |12| . 
For P within the disk around the origin we refer to [3~8] Section 4]. 

Note that p(°°) and S have the same jumps on J \ {0}, and that P and S have the same 
jumps on the lens X within the disks. This implies that R is analytic on the entire plane, except 
for jumps on the reduced system of contours Hr, as shown in Figure H3 cf. 12 , and except for a 
possible isolated singularity at the origin. Yet, as in |27| 138 j. it follows easily from the behavior 
of S and P near the origin, given by ()3,12|) and (|3.13|) , and by (|5.4|) and 1)5.5)1 , respectively, that 
the isolated singularity of R at the origin is removable. Therefore R is analytic on C \ T,r. 

Recall that the matrix valued functions S and p(°°) are normalized at infinity. Since 
detP^ 00 ) = 1, this implies, by (|6.2[) . that also R is normalized at infinity. 

Let vr be the jump matrix for R. It can be calculated explicitly for each component of E#. 
However, all that we require are the following estimates, cf. 

||w R (z)|| = 1 + 0(e" cn|z| ), as n -> oo, z £ T, R \ circles, 
1 1 1 1 = I + 0(l/n K ), as n — > oo, z G circles, 

for some c > and < k < 1, and where || • || is any matrix norm. We note that the extra factor 
|x| 2q , which we will meet in vr, does not cause any difficulties to obtain this behavior. These 
estimates then imply that vr is uniformly close to the identity matrix as n — > oo, and, since R 
is normalized at infinity, we then find uniformly for z E R \ T,r, 



R(z) = I + 0(l/n K ), asn^oo. (6.3) 

d_ 

dz J 



So, R is uniformly bounded as n — > oo. We also have that 4zR(z) is uniformly bounded as 



n — » oo. Another useful property is detP = 1, which follows from 1)6.2)1 and the fact that 
S p(°°) and P all have determinant 1. 
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7 Proof of Theorem 11.11 

We now have all the ingredients necessary to prove Theorem ll.il We point out that the general 
scheme of this proof is the same as the proof of |28[ Theorem 1.1(c)]. We replace in the kernel 
K n , given by (|1,3|) . the orthonormal polynomials p n -i,n an d p n .n, together with their leading 
coefficients 7 n -i,n and 7 n ,n, by the appropriate entries of Y, given by (|2.4|) . and find 

;Y n (x)Y 21 (y) -Y 21 (x)Y u (y) 



K n (x,y) = - — \/w n (x)y'w n (y) . (7.1) 

lui x — y 

This means that the kernel K n can be expressed in terms of the first column of Y. Hence, we 
want to know the asymptotic behavior of Y near the origin. This will be determined in the 
following lemma. 

Lemma 7.1 For x S (0,5), 

r„(,A _ fT-,^) , (7 . 2) 



wzi/i M(z) given by 

M(z) = R(z)E(z)e n ^°^e-^^= Q , (7.3) 

where R is the result of the transformations Y^T^S^Rof the RH problem, and the 
matrix valued function E is given by (|5.27j) - (|5.3()|) . The matrix valued function M is analytic 
in Us with M(z) and 4-M(z) uniformly bounded for z G Us as n — > oo. Furthermore, 

detM(*) = l. (7.4) 

Proof. We use the series of transformations Y — ► T —* S —* i? and unfold them for z inside 
the disk t/5 and in the right upper part of the lens, so that z 6 /~ 1 (i'). Since w(z) = z 2a and 
W(ar) = z a e- TiQ for our choice of z, see (jSS| and (jBTSjl . we have by tiTTty . (pT7]l . (ET251 and 



z -2a e -2n</»( 2 ) \J 

We then get for the first column of Y , 



(7.5) 



SiS) = z " aen(s{2) ~ 0w_ ^ e ^ a3jR (' z ) £; «( z )^( n /( z )) e ™ <T3 (1) • ( 7 - 6 ) 



Since z is in the right upper part of the lens and inside the disk Us, we have < arg nf(z) < 7r/4, 
cf. Figure and we thus use (|5.23[) to evaluate ^f a (nf(z)). Using the formulas 9.1.3 and 9.1.4 
of which connect the Hankel functions with the usual J-Bessel functions, we find 

/1\ ™ f(nf(z)) 1/2 J n ,i(nf(z))\ 

* (n/(z))e^ (I) = e-x^ ' a+ * 7 . (7.7) 

By (|5~26|) and (|7~3|) we have R(z)E n (z) = M(z). Inserting this and <|T7J) into (|7~6|) we get 

An(,)\ = e -^^, e ^ (<) f (n/W)VV a+| (n/(z))\ 
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Letting z — » x £ (0,5), and noting that 



x e 



°/- P n{g+(x)-4>+(x)-^) _ x ~a e \n(g+(x)+g-{x)- 



a \nV{x) 



X e2 



w n (x)' 



-1/2 



(7.9) 



which follows from the fact that 2<p + (x) = g+(x) — g-(x), see Sectional and from (|1.2j) and 
(t2~TUl) . we obtain ([7^1 . 

The matrix valued function M is analytic in the disk Us since both R and £7 are analytic in 
this disk. So, we may write M(x) instead of M + (x) in l|7.2jl . 

We recall that R(z), jzR(z), E(z) and -^E(z) are uniformly bounded for z £ f/j as n -> oo, 
see Section El and Section HO If we also use that |e n ^+(°)| = 1, which follows from the fact that 
0+ is purely imaginary on J, we have from Q7.3|) that M(z) and 4-M(z) are uniformly bounded 
for z £ Us as n — ► oo. 

Since M is a product of five matrices all with determinant one, ()7.4|) is true. □ 
Lemma 7.2 Let u £ (0, oo), u n = n ^ ) a^d = nf(u n ). Then 



II,, = TTU + O [ ^— ) , 



as n —* oo, 



(7.10) 



as n — ► oo, 



(7.11) 



^ Q _l(^n) = ^ Q _l(vrn) +0 

2 2 \ n 



as n — > oo, 



(7.12) 



where the error terms hold uniformly for u in bounded subsets of (0, oo). 
Proof. Since, see (|5.8|) 



f(x) = 7T^(0)X + 0(X 2 ), 



as x — > 0, 



we have, uniformly for it in bounded subsets of (0, oo), 



/ 



— 1 



u 



TT- + 



n 



u 



as n — > oo, 



which proves (|7.10j) . 

We note JTJ formula 9.1.10] that J Q ,i(z) = z a+ ^H(z), with H an entire function. It then 
follows from (|7.10|) that, as n — > oo, uniformly for u in bounded subsets of (0, oo), 



(^) Q+ 2 + O 



V+i 
n 



H(iru) + O — 
n 



so that equation Q7.11JI is proved. Similarly, we can prove (|7.12|) . 
We are now able to prove Theorem ll.il 



□ 
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Proof of Theorem II. 1L Let u, v € (0, oo) and define 



u 



u„ 



nV(0)' 



v„ 



nf(u n ), v n = nf(v n ). 



We put 

K n (u,v) 



1 



-K n (u n ,v n ). 



From (|7,1[) and (|7.2|) we then have 

1 



K n (u,v) 



2iri(u — v) 



det 



Ttz i r~ 

)Yii(u n 



e 2 



Un ) 



2(u - v) 



det 



M(u n ) 



'-1/2 T ( ~ \ 
Un J a+ i(U n ) 

Un 2 J a _i(u n ) 




i \ V r , 



The matrix in the determinant can be written as 

uj 2 J a+ ±(u n ) vl /2 J a+ i(v n ) 

Un ^a-liUn) V n /2 J a _i(v n ) 



M(v n ) 



+ M(v n )- 1 (M(u n )-M(v n )) 



Un 2 J a+ i(u n ) 0' 
Un 2 J a _l{u n ) 



(7.13) 



We will now determine the asymptotics of the second term in (j7.13j) . Since detM(v n ) = 1 
and since ikf(z) is uniformly bounded for z £ see Lemma 17. l| the entries of A£(v n ) are 
uniformly bounded. By Lemma 17. II we also have that -^M(z) is uniformly bounded for z E Ug, 
so that from the mean value theorem M(u n ) — M(v n ) = O (^-^)- From Lemma 17.21 it follows 

that u n J , i (u n ) = 0(u a+1 ) and ul/ 2 J i(u n ) = 0(u a ) uniformly for u in bounded subsets 

' 2 2 

of (0, oo) as n — > oo. Hence we have, uniformly for u, v in bounded subsets of (0, oo), 



M{v n )-\M{u n )-M{v n )) 



u n /2 J a _du n ) \0 0, 



1 /2 

Inserting this into (|7.13() . using the fact that detM(v n ) = 1, and that v n J a ±i(v n ) = 0{v a ) as 



n — > oo, we then find uniformly for u, v in bounded subsets of (0, oo), 

'"M 2 J a+ 1 («n) + O {^U a ) Vn /2 J a+ l («„) 1 

.1 
Ur. 



K n (u,v) 



1 



det 



2(u-v) \^j a _ k ^ n ) + 0(^u a ) vl /2 J a _i_(v n ) 



1 



Un 

.(u n ) Ih/ 2 J i{v n ) 
2 2 



2(u - v) 



det 



-1/2 T 



+ o 



r? 



(7.14) 



We note, from Lemma 17.21 that we can replace in the determinant, u n by ttu and v n by irv. 
We then make an error which can be estimated by Lemma 17.21 However, since this estimate is 
not uniform for u — v close to zero, we have to be more careful. We insert a factor u~ a in the 
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first column of the determinant, and a factor v a in the second. Then we subtract the second 
column from the first to obtain 



K n (u,v) 



U a Un 2 J a+ l (U n ) ~ V a vl/ 2 J a+ l(v n ) V a vl /2 J a+ l(v n ) 
2{ U - V) ^ \u- a vt/ 2 J a _i (U n ) - V- a V l n /2 J a _i_(v n ) V- a Vn /2 J a _i(d n ) 



det 



+ 



n 



(7.15) 



One can check, using (|7.1t)j) . (|7.11jl and the facts that J' a+ i(x n ) = J' a+ i ( 7rx ) + 0{ 
j-x n = vr + 0(f), where we have put x n = nf{-§ m ), that 



and 



dx 



ntp(0) - 



d r 

dx 



- a xH 2 J a+k {x n ) - x~^x)^J a+l {vx)\ =0Q, 



as n — > oo, 



uniformly for x in bounded subsets of (0, oo). It then follows that the (l,l)-entry in the deter- 
minant of (|7.15j) is equal to 



u 



- a {nu) l l 2 J a+ ,_ (vrn) - v^iTTv) 1 / 2 J a+ i (nv) + O ( 
Similarly, we have from 



u — V 



n 



d_ 

dx 



oil 

n 



as n — > oo, 



u — v 



that the (2,l)-entry in the determinant of (|7.15|) is equal to 

u- a (TTu) 1/2 J i(vru) -v- a {itv) l ' 2 J n _ ihrv) + ( 
2 a V 

From Lemma l7.2l it also follows that v n J a ±l (vn) = ( 7 tv)J a ±i(Trv)+0(l/n). Therefore, uniformly 
for u, v in bounded subsets of (0, oo), 

K n (u,v) 



x det 



2(n - v) 

'u-ot-Ku) 1 ' 2 J a+ i (vrn) - v-^irv) 1 ' 2 J a+| (ttv) + O (^) iT^™) 1 /^! M + O 



+ o 



J£(u,u) + 



+ o 



n 



2(n - v) 



det 



it-H^^i (ttu) - ^(vrw) 1 ^ j Q _i (ttu) O (I) 



Since z a+ 2 J a± i (z) is an entire function we have by the mean value theorem that 

u- a (TTu)^ 2 J a± i (vrn) - ^-"(Tru) 1 ^ j a± i (ttu) 
u — V 



(7.16) 
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is bounded for u,v in bounded subsets of (0, oo). From (j7,16|) we then have 
K n (u,v) = S° a (u,v) + 

uniformly for u,v in bounded subsets of (0, oo), which completes the proof of Theorem ll.il □ 
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